Cooling a magnetic resonance force microscope via the dynamical back-action of 
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We analyze the back-action influence of nuclear spins on the motion of the cantilever of a magnetic 
force resonance microscope. We calculate the contribution of nuclear spins to the damping and 
frequency shift of the cantilever. We show that, at the Rabi frequency, the energy exchange between 
the cantilever and the spin system cools or heats the cantilever depending on the sign of the high- 
frequency detuning. We also show that the spin noise leads to a significant damping of the cantilever 
motion. 
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I. INTRODUCTION 

Magnetic resonance force microscopy (MRFM) 
is a powerful technique for visualizing subsurface 
structures^i^i^i^ with three dimensional spatial reso- 
lution of the order of 10 nanometers or less^, which 
is more than two orders of magnitude better than the 
resolution of conventional high-field magnetic resonance 
imaging (MRI). The main part of a MRFM device is 
a nanomechanical resonator or cantilever (see Fig. [T] 
below), having a fundamental frequency in the range 
of several kHz. By using MRFM, a significant break- 
through in magnetic resonance detection sensitivity was 
achieved, resulting in single-electron spin detection^ with 
a spatial resolution ^25 nm and substantial progress 
in nuclear spin detectiort^'^ i-'^^'^^i^^ . MRFM has also 
been proposed as a qubit readout device for spin-based 
quantum computers^^ii^. 

MRFM was initially proposed as a possible means 
to improve the detection sensitivity to the single spin 
leveli. Since then, progress in MRFM and related 
technologies has attracted broad interest, especially the 
questions of squeezed states of the cantilever and the 
collapse of its wave function when both, the spin to 
be measured and the cantilever, are treated quantum 
mechanicallyi^ii^iiiiiSiii^. However, the ultimate goal to 
detect a single nuclear spin with MRFM is still a chal- 
lenge. Rough estimates indicate that in order to reach 
this goal the effective temperature of the MRFM can- 
tilever should be reduced to about 0.1 /iK, which corre- 
sponds to ~2.5 kHz of the fundamental frequency of the 
MRFM cantilever. 

Numerous experiments on cooling micro-mechanical 
resonators via their coupling with different external sys- 
tems have recently been reported (see, e.g., Ref. [20l . 
[2lll2^l23ll2^ . Experimental results show that a micro- 
mechanical resonator can be cooled down to an effec- 
tive temperature on the order of 0.1 K(Ref. [20l ) or 5 
mK (Ref. [23l ). However, in order to drive the micro- 



resonator to the quantum regime, more effective cool- 
ing methods are needed. A promising way would be to 
cool the micro-resonator by coupling it to a solid-state 
quantum electronic circuit. In principle, the effective 
electronic cooling of the micro-resonator can be achieved 
by several means, including coupling it to an another 
resonator.25i2S, to a transmission line resonator^L^^, to a 
quantum dot^, to an electronic spin'^*', or to supercon- 
ducting dcvicc a^^'^^-^^-^"^-'^^-^^'^'''^^'"^^'"^^''^^'"^"^'"^^ 

In particular (as shown, e.g., in Refs. [35ll37ll39ll40ll4ll . 
I42i i43.44i ) , an electric resonator circuit weakly coupled to 
a two-level system (superconducting flux qubit) can be 
cooled by its quantum-dynamical back-action. 

In this paper, we investigate the cooling of a MRFM 
cantilever via its coupling to nuclear spins. We show that 
the back-action of the spin system modifies the equation 
of motion of the cantilever, providing additional damping 
and a frequency shift which depends on the properties of 
the spin system (decoherence rates, damping rates, etc). 
We investigate the operation modes of the MRFM where 
the damping is positive and results in a substantial de- 
crease of the effective quality factor of the MRFM can- 
tilever, and thus a significant cooling of the cantilever 
motion. 

This paper is organized as follows. In Section [H] we 
describe the interaction of a MRFM cantilever with nu- 
clear spins. We obtain the cantilever equation of mo- 
tion modified by the back-action of the spin system. The 
modification appears as an additional contribution to the 
damping and frequency shift of the cantilever, in terms 
of the magnetic spin susceptibility. 

In Section IIIII we obtain the explicit expression for the 
low-frequency spin susceptibility. The most important 
result of this section is that the longitudinal magneti- 
zation of a sample has a clear resonance at its Rabi fre- 
quency. In some sense, this is the low-frequency analog of 
the conventional high-frequency NMR for the transverse 
magnetization. 

Section IIVI is devoted to a detailed study of the influ- 



ence of the spin system on the damping and the frequency 
shift of the cantilever. In the first part of this section we 
consider a sample with a relative short spin-lattice re- 
laxation time Ti. In this case, we show that the quality 
factor of the cantilever changes depending on the sign 
of the high-frequency detuning. For positive detuning 
(the microwave frequency is above the nuclear resonance 
frequency), the contribution of the spin system to the 
cantilever damping is negative: heating the cantilever by 
absorbing Rabi photons from the spin system. If the 
detuning is negative (the microwave frequency is below 
the nuclear resonance frequency) , the contribution of the 
spin system to the cantilever damping is positive: cool- 
ing the cantilever by giving up Rabi photons to the spin 
system. In the second part of Section IIVI we consider 
the influence of the spin noise on the cantilever motion. 
Independently of the particular values of the parameters 
which characterize the spin system, we show that its in- 
fluence on the cantilever damping is always positive, i.e., 
the spin noise always leads to a decrease of the quality 
factor of the cantilever. 
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II. INTERACTION OF NUCLEAR SPINS WITH 
THE MRFM CANTILEVER 

A schematic diagram of the system studied here is 
shown in Fig. [TJ A spherical ferromagnetic particle with 
magnetic moment mp is attached to the cantilever tip. 
A small paramagnetic cluster with magnetic moment /x, 
which must be detected, is placed on the surface of a 
non-magnetic sample beneath the tip of the cantilever. 
The whole system is placed in a permanent high mag- 
netic field, Bq, oriented in the positive z-direction. The 
transverse magnetic field Bi{t), which excites the NMR 
in the sample, is applied to the paramagnetic cluster. In 
addition to Bq, the magnetic moment /i experiences the 
inhomogeneous field Bp{z) from the ferromagnetic tip. 
We assume that the field Bp{z) is also oriented in the 
positive z-direction and is given by the dipole formula: 



Bp{z) 



fiQ 2mp 
Att {d + zY 



(1) 



where d is the equilibrium distance between the cantilever 
and the sample surface, and z is the amplitude of the 
cantilever oscillations. Below we assume z d, hence 



FIG. 1: Schematic diagram of the setup of the system under 
consideration. Here, -Bo is the uniform permanent magnetic 
field, Bi{t) is the rotating rf magnetic field, F{t) is an external 
force acting on the cantilever in the ^-direction, mp is the 
magnetic moment of the ferromagnetic particle attached to 
a free end of the cantilever, fj, is the magnetic moment of 
the nuclear spin, located in xy- sample plane, and d is the 
equilibrium distance between the center of the ferromagnetic 
particle and a sample plane. 



The interaction of a particle, having a magnetic mo- 
ment /z, with the cantilever is given by the following 
Hamiltonian 



H = He + Hq - 



- HcB + Hs + H^g'^ - 



H. 



SB - 



-at-B), (3) 



where He is the Hamiltonian of the cantilever 



He = 



2™ 



(4) 



Hs is the Hamiltonian of a nuclear spin interacting with 
a one-mode high-frequency field 



Hs 
where 



~^'^z ^ — T^'^xia + a) + hujni^a^a , (5) 



Bp{z)^Bp{0) ( 



where 



BpiO) 



2^omp 



(2) 



Hence, the cluster we investigate is under the polarizing 
field: Bq + Bp{0). 



= -f{Bo + Bp{0)) 



7 is the nuclear gyromagnetic ratio, Wmw is the frequency 
of the microwave field. The quantity hui, the interaction 
energy between spin and the microwave field, is propor- 
tional to its amplitude Bi (see below). The magnetic 
moment 7? of a spin-1/2 particle is expressed in terms of 
the Pauli spin matrix vector: 'jt — —h'flj /2. 



The Hamiltonians H^^ and H^' represent the ther- 
mal baths for the cantilever and spin, respectively, while 
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HcB and Hs b represent their interactions with their cor- 
responding baths. 



We will not specify here the bath Hamiltonians i/, 



(b) 
C ' 



Hg'\ and their interactions Hcb, Hsb- We describe the 
influence of He, Hcb on the motion of the cantilever by 
introducing the damping rate 7c and the external noise 
r,{t). 

We consider the cantilever tip as an oscillator with 
effective mass m, and effective spring constant kg, subject 
to an external force _Fo cos wt and force fluctuations with 
a spectral density 



muJi 

~Q'c 



■ coth I ■ 



\2kT 



(6) 



where Qc is the quality factor of the bare cantilever. 

The equation of motion for the cantilever, interacting 
with N spin-1/2 particles, then reads: 



Mz dBp 

ni dz 



+ focosi^t + r]{t), (7) 



external force f{t). It has the magnetic susceptibility 
Xzzi'-^) of the spin system as its Fourier transform: 



Sfiti) 



duj 
2^ 



c^p[-iuj{t - ti)]xzz{(^) , (10) 



From Eqs. ([9]) and ([8]) we obtain the equation of motion 
of the cantilever: 



+ X^m / dti 



2- _ ^/_(o)(,^\ 4A/^(o) 



z + 7ci + uj^z = \{af {t)) - - (af {t)) z 



Sfih 



z{ti) + focosi^t + Tj{t) , (11) 



Let us now analyze this equation in detail. The first two 
terms in the right-hand side of Eq. pT|) contribute, re- 
spectively, to the amplitude and the frequency shifts due 
to the spins. The modulation of either of these terms is 
usually employed in MRFM experiments. In particular, 
the second term of the right-hand side of Eq. pT|) can be 
easily converted to the frequently used expression for the 
frequency shiftiS: 



where Adz is the quantum operator of the longitudinal 
magnetization, 

Mz = —(Jz , 

ujc = (ks/m)'^/'^, -fc = LOc/Qc, fo = Fq/tu, and ri{t) de- 
scribes the fluctuations of the acceleration, and has a 
spectral density S',,(w) = Spiuj)/!!!?. 
By using Eq. ^ we obtain: 



Mz dBp 
m dz 



A(crz) - ^{zaz) 



(8) 



where 



A = 



2md 



is the coupling strength between the spin and the can- 
tilever. The angular brackets in ^ and ^ denote the 
average over the two free baths variables. 

The quantity (ctz) in ([5]) is a functional of the can- 
tilever position z(t), where its first order (« z/d) contri- 
bution to the Hamiltonian ([3]) is {Xmz + f)<Tz, where we 
introduce f(t): a small external force that is required for 
calculating the magnetic susceptibility^. Hence, to first 
order in A, we obtain 



where (cr^'' (t)) is described by the evolution of the spin 

system {Hs + H^'^ +Hsb) uncoupled from the cantilever. 

The functional derivative 6{az{t)) /Sf{ti) in Eq. © is 
the response of the spin system to the weak low-frequency 



.(0) 



(t)) + Xm 



S(azit)) 



(9) 



A/ = mJ,{d^B,/dz^)/2k, , 

where iriz is the magnetic moment of the sample. These 
two terms describe the direct influence of the spin on 
cantilever motion. 

The third term in the right-hand side of Eq. (fTTj) de- 
scribes the additional back-action of the spin on the can- 
tilever motion, which is due to the modification of the 
spin dynamics by the cantilever. This term, which is the 
main subject of our study here, gives rise to an addi- 
tional frequency shift and an additional damping of the 
cantilever. 

Assuming that the steady-state value of the quantity 
{o''"z\t)) is independent of time, we convert Eq. (Ilip to 
Fourier components of z(i) [using z(i) = / e~*"*2;(w)(ia;)]: 



{ 



^c^--^ + ^(4"Vaw.(-) 

- i [w7c -t- X'^mx'^.iuj)] |z(a)) = y(5(w - i^) + r]{uj) , 

(12) 

where we introduce the real and imaginary parts of the 
spin susceptibility 

Xzzi^^) = X'zzi^) + ^Xzzi^) ■ 

Analyzing the third and the fourth terms of the rhs 
of Eq. we see that the influence of the spins on the 
cantilever produces the frequency shift 



dujc ^ 2uJc 



(13) 



where the flrst term in the rhs of Eq. represents the 
direct contribution of the spins to the frequency shift. 
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while the second term in the rhs of Eq. (|T3)) is the ad- 
ditional contribution to the frequency shift that results 
from the indirect influence of the back-action of the spins 
on the cantilever motion. 

From the fifth term in the rhs of Eq. (fT2)) . 
— i [a;7c + A^r7ix"^(cj)] z(ct;), we can write the total 
damping of the cantilever: 



Ttotal Tc ~^ Tsi 



(14) 



where 7spin is the frequency-dependent contribution of 
the back-action of the spins to the damping of the can- 
tilever: 



(15) 



From this expression we obtain a spin back-action- 
induced modification of the cantilever quality factor: 



1 _ 1 



(16) 



III. THE LOW-FREQUENCY MAGNETIC 
SUSCEPTIBILITY OF IRRADIATED SPINS 



As was shown in Ref. l43l this system can be described, 
within the RWA approach, by the rate equations for the 
elements of the reduced density matrix, which describe 
the transition between Rabi levels [the levels £'^(7Vph) in 
Eq. ini), with the same iVph]. 



dp 
'dt 



= -Aip + Bp+ + (r_)cos26' 



(19) 



dp+ 
dt 



-iniip^ + Bp- A2P+ + iT^)sm20 , (20) 



where 



~dr 



A. 



— cos^ 29 + sin^ 261 



— sin^ 29 + cos^ 261 
J-i 



(21) 

(22) 
(23) 



The interaction of a two level system with an exter- 
nal electromagnetic field tuned near the resonance of 
this two-level system can be described using the dressed- 
state approach^, which recently was successfully applied 
to investigate the interaction between solid state super- 
conducting qubits and an external radiation source^^ii^. 
Some results obtained in Ref. within the dressed-state 
approach will be applied here for the investigation of the 
interaction of a MRFM cantilever with irradiated nuclear 
spins. 

The energy levels of a spin-1/2 interacting with a high- 
frequency Wniw field has a term proportional^ to the 
number A^ph of photons, with an additional splitting of 
each photon state by the Rabi energy Mlfi/2: 



where riu is the Rabi frequency 



(18) 



with 



where (A^ph) is the average number of high frequency 
photons^^, and 



- too 



is the high-frequency detuning. For definitiveness, here 
we assume S > 0. The frequency f2i is directly related to 
the amplitude of microwave field: f2i = jBi. 



B 



sin 29 cos 26* 



(24) 



where F;^ is the dephasing rate of a spin, which can be 
expressed in terms of the spin-spin relaxation time T2, 
F;^ = 1/T2. Here, Ti is the spin-lattice relaxation time, 
which is related to up F| and down F| transition rates 
between spin levels 

T,-' = Ft + Fi ; F_ = Ft - F^ . 

The angle 9 is defined by tan26' = —ili/S, where < 
26* < TT, so that 



cos 261 = -6/nR , 



(17) and 



cos9 = —=[ 1 - — 

V2\ 



1/2 



sin.„(l 



1/2 



For equilibrium conditions, the relaxation F| and ex- 
citation F| rates are related by the detailed balance law: 



Tt = Tx exp 



From Eq. 



we obtain 



TiF- 



tanh 



hiuo 



^0 
2fcsT 



(25) 



(26) 



The quantity p in Eq. is defined as the difference of 
the populations between the higher and the lower Rabi 
levels. 
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The steady-state solution = 
Eqs. (Ull), ® , and ^ is as follows: 



dt 



(r: 



(r_) COS 261 , 



r2 



■(r_)sin2e' 



0) for 



(27) 



(28) 




P 



(0) 



.^^fl (0) 



(29) 



It is interesting to note that under high-frequency ir- 
radiation, the population of the Rabi levels becomes in- 
verted. This is seen from Eq. (j27p . where the quantity 
p^"^ is positive, since for 5 > we have cos 20 = —6/ flu < 
0, and always r_ < 0. 

In addition, as S tends to zero, p*^") 0, which causes 
the equalization of the population of the two levels when 
the high frequency irradiation is in exact resonance with 
the NMR frequency jBq. 

The quantity {az{t)) which is, by definition, the longi- 
tudinal magnetization of a sample with TV spin-1/2 par- 
ticles (see Appendix A) can be expressed^ in terms of 
the matrix elements p and p+ : 



{(Tzit)) = p{t)cos29 + p+{t)sin2e 



(30) 



Therefore, from the definition of Xzz{^^) in Eq. PH)) we 
obtain 

Xzz(w) = Xp{^^) cos2e + Xp+{uj) sin2e , (31) 

where Xpi^) ^nd Xp+i'^) ^tre the spectral components of 
the response of p{t) and p+{t) to a weak external force, 
and are defined similarly to Eq. pH]) . 

The susceptibilities Xpi'^) ^nd Xp+i^) can be readily 
found by investigating the response of the reduced den- 
sity matrix in Eqs. (|19p. ^U\i . and (PT|) to a weak external 
perturbation [see Eqs. (|B8p . (|B9|) . and (jBlOp in Appendix 
B]. Hence, the expression for Xzz{^) becomes: 



Xzzi^^) 



6 r^i fiij 



pf (-icj + 2r^) 



(32) 



where D{Ld) is given in Eq. (jBlip . 

The quantity p^"* in Eq. ([32|) is the steady-state value 
which can be written as follows: 



(0) 



1 



nRl + {5T2Y+T2T^nl 



tanh 



hijjQ 
2kBT 



(33) 



and can be expressed in terms of the stationary mag- 
netization (jAlOp : 



Nh-y (0) 
^^P+ = 



M 



(st) 



1 + (<5r2)2 nn 



(34) 




B,cosco2t 



FIG. 2: Resonance of the longitudinal magnetization Mz 
at the Rabi frequency Qr. (a) For conventional NMR, 
the longitudinal magnetization of a sample polarized in the 
field -Bo and subject to a circularly polarized excitation field 
Bi cos LUit, where wi ~ wo, exhibits a precession around the z- 
axis with a time-independent steady-state ^-projection of the 
longitudinal magnetization M^. (b) A second low-frequency 
excitation B2 cosuj2t, where W2 ~ ^r, applied along the z- 
axis, produces resonant oscillations of Mz near the Rabi fre- 
quency Qr. 



A. Resonance of the longitudinal magnetization at 
the Rabi frequency 

Historically, the detection of NMR is based on the 
Faraday law of induction^". That is why most of the 
measurement schemes in NMR are based on the detection 
of the transverse magnetization, which oscillates with a 
relatively high frequency. 

The detection of the longitudinal magnetization is less 
common because it requires measurements in a low fre- 
quency range with a low signal-to-noise ratio. How- 
ever, this drawback can be circumvented by some tech- 
niques, such as the pre-polarization of a sample in high 
field^, or the use of superconducting quantum interfer- 
ence devices (SQUIDs)'"'^, which allow to obtain, in the 
micro- and nano-Tesla range, a resolution which is be- 
yond what is usually achieved in conventional high-field 

j^jjjp 53 .54.55.56.57 

In MRFM there is no choice other than to measure the 
longitudinal component of the nuclear polarization, since 
the resonance frequencies of MRFM cantilevers are well 
below those corresponding to the frequencies of NMR 
transitions. From this point of view, it is interesting to 
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note that the longitudinal magnetization of a sample 
placed in a high frequency resonant radiation field shows 
a clear resonance at the Rabi frequency if the sample is 
subject to an additional low-frequency excitation directed 
along the z-axis with energy (see Fig. [2]) 

Hlf = ^7i?2 cos ojt . 

As shown in Refs. and m, this effect is a general 
feature of any two-level dissipative system. In this case, 
the low-frequency evolution of the longitudinal magneti- 
zation can be expressed in terms of the spin susceptibility 



Xzzi^) cosa)t+x'4(w) sin ujt 
(35) 



where * is given by Eq. (|AlOl) . 

With the aid of Eq. 
following form: 



Mzit) 
where 



M 



(st) 



1 



Bi 



we rewrite Eq. (|35p in the 

[Xzz (^) COS ujt + Xzz {^) sin uit] 



DiLu)Tl{l + {T2Sy) ' 



(36) 
(37) 



Bi is the amplitude of the /ii(?/i-frequency resonance ex- 
citation, and B2 is the amplitude of the Zow-frequency 
signal which excites the Rabi oscillations of the longitu- 
dinal magnetization M^- 

The resonance at the Rabi frequency is clearly seen 
in Fig. [31 where we plot (for two values of the spin-spin re- 
laxation time T2) the dissipative part of the spin suscepti- 
bility p7p as a function of the low- frequency uj [uj fa ^r)- 

It should be remembered that our linear approximation 
is valid within the range 



(38) 



from where we obtain the range of the amplitudes of the 
low- frequency signal B2 , where the expression (|36)) is con- 
sistent with the condition (l38ll: 



B2 < Bi/x 



If 

zz. max 



(a.) 



Hence, near the Rabi resonance (i.e., to ~ ^r), a sig- 
nificant modulation of the longitudinal magnetization Mz 
can be induced by a low frequency drive. 

The resonance of the longitudinal magnetization de- 
scribed above is, in some sense, a low-frequency analog 
of conventional NMR, where the resonance of the trans- 
verse nuclear magnetization is being measured at a high 
resonance frequency. 



-50- 



3, 

N 
N 



-100 



-150 
0.' 



^^^^^ 




T, = io-sN\ 


It 




- 




B„ = 3T 




= 0.001 T 




= 0.1 s 




5/cOp = lO""* 


= 2x1 0"'' s I / 









1.1 



FIG. 3: (Color online) The dependence of dissipative part 
^"^(aj) of spin susceptibility for low frequencies lo ~ f2_R. 



IV. INFLUENCE OF SPINS ON THE 
FREQUENCY SHIFT AND THE DAMPING OF 
THE CANTILEVER 

A. Influence of the driven equilibrium longitudinal 
magnetization on the damping and frequency shift 
of the cantilever. 

Here we consider the situation where, before the mea- 
surements start, the spins under the application of a po- 
larizing external magnetic field reach thermal equilibrium 
with their environment. In other words, in this case, the 
spin-lattice relaxation time T\ is sufficiently short (for 
example, in the millisecond range) to ensure the applica- 
tion of conventional NMR measurement protocols. For 
this case, the corresponding susceptibilities are given by 
the expressions (IB12|) . (|B13I) . and (|B14p . 

In order to analyze the damping and the frequency 
shift of the cantilever in Eq. we explicitly write 

down the real and imaginary parts of Xzz{^) from 
Eq. — 



(39) 



x'zzH 

where 



D,{Lo)=Re[DiLo)] = —[AiTinl + Tl~Lo^{l + 2TiT^)] , 

(41) 

D2iu;) ^ MDiu;)] = (^u;^ - - _ . (42) 

It is worth noting that at the exact resonance {S = 0), the 
back-action influence of the spins on the cantilever van- 
ishes, since the corresponding susceptibilities and 
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(|40| are equal to zero at this point. This is a conse- 
quence of the equahzation of the population of Rabi lev- 
els at the point of exact resonance (p^"-* tends to zero as 
S approaches zero). This produces the vanishing of the 
energy flow between the spins and the cantilever at the 
Rabi frequency. 

Another point is that the sign of the susceptibility 
Xzzi^) Eq. (|40p is opposite to that of 5, for any value 
of its parameters. This follows from the fact that the nu- 
merator ujDi{uj) + 2r^D2{uj) in Eq. (|40)) is always nega- 
tive, and the quantity p^"* is also negative. Therefore, for 
S > (when the higher Rabi level is more populated than 
the lower) the contribution of the spins to the damping of 
the cantilever is negative [see Eq. fT^ ]. In this case, the 
cantilever is being heated by absorbing the Rabi photons 
emitted by the spins. 

In the opposite case, when S < (i.e., the higher Rabi 
level is less populated than the lower) the contribution 
of the spins to the damping of the cantilever is positive, 
therefore cooling the cantilever which gives up Rabi pho- 
tons to the .spin system. 

It is worthwhile to consider the dependence of the dis- 
sipative part of susceptibility Xzzi^) (Eq- (HO])) on the 
frequency to. It has two peaks. A lorentzian peak is in 
the vicinity of the Rabi frequency flu, as it is evident 
from ([1^ . The condition for this is the relative large 
high-frequency detuning (S 3> F^, fii). The approximate 
expression for the peak value at the Rabi resonance is as 
follows: 



peak 



(0) 



2hfl 



(43) 



The other peak is related to the spin-lattice relaxation 
time Ti of the longitudinal magnetization. It lies at much 
lower frequencies (w <C ^Ir). If we assume ui ^ flujT^p; 
S > r^,Ti,f]i; r^Ti > 1, we obtain Di{lj) « 6'^/Ti, 
D2{io) « —S'^oj, which yields the following expressions 
for the susceptibility: 



x'zA^) 



.(0) 



5 2VLiTi 



1 



\5\ hT^ 1 + 



(0) S 2fliTi LoTi 
\5\ hT^ 



(44) 



(45) 



with the maximum of Xzzi^) being at Wmax = 
These expressions are analogous to the Debye formu- 
lae for the low-frequency dispersion of the dielectric con- 
stant. In the context of NMR, this Debye- like behavior is 
known for the response of the transverse magnetization 
in a weak polarizing field^'^. 

To estimate the magnitude of this effect, we take as 
a guide the experimentally accessible parameters from 
Ref . IT^. From Fig. 2b of that paper we take the following 
values of the magnetic field Bp{0) on the tip, Bp{Q) = 
100 mT, and a corresponding value of d, the distance 
between the tip and the sample: d = 60 nm. From the 
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FIG. 4: The imaginary part of spin susceptibility as func- 
tion of normalized frequency uj/Qr in the vicinity of Rabi 
frequency, for T = 0.6 K, N=10^ 



spring constant kc — 6 x 10^ 
frequency of the cantilever, fc 
cantilever mass 



N/m and the resonance 
= 3 kHz, we estimate the 



TO 



kjuj^^ = 1.6 X 10"'-^ kg. 



For the element ^^F, which was the subject of study in 
Ref. [13, the nuclear gyromagnetic ratio 7 = 27r x 40 
MHz/T. This allows to estimate the coupling factor A: 



A = 



2md 



4.12 X 10" 



Here N is the number of spins in the resonant slice. Thus, 
for the factor X^m/h we obtain the estimate 



(2.6 X 10* X TV 



10- 



In addition, we take Ti = lO""^ s, T2 = 1/F^ 
T = 0.6 K, Bo = 3 T, Bi = Qi/j = IQ-^ T. 

The contribution of the spins to the cantilever damp- 
ing 7spin is proportional to the imaginary part of the 
spin system susceptibility x'Li^) i^^^ Eq. p^ ]. The 
susceptibility Xzzi'^) has a clear resonance at the Rabi 
frequencies, which is shown in Fig. ID The second Debye- 
like peak is shown in Fig. [S] It is worthwhile to note 
that for the parameters we used here the peak value in 
Fig. [5] is much higher than the Rabi peak in Fig. |4l 

These resonances modify within the corresponding fre- 
quency range the bare cantilever quality factor Qc in 
Eq. (fTB)) , as shown in Fig. [H] and Fig. [T] 

The contribution to the cantilever frequency shifl is 
given by the real part of the spin system susceptibil- 
ity x'zzi^) Eq. ([55)1 . The associated frequency shift, 
A^TOXzz('-^)/2wc, near the Rabi resonance is shown in 
Fig. [H and near the Debye-like peak in Fig. [9] 

For the parameters used here, the Rabi frequency 
ilu w 12 MHz and the Debye-like peak in Fig. [5] is lo- 
cated near 160 Hz. Between these two peaks, the quality 
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FIG. 5: The imaginary part of spin susceptibility as a func- 
tion of the normahzed frequency u/Qr, in the vicinity of the 
Debye-like peak, for T = 0.6 K, and N=10^ 




FIG. 6: The modified quality factor Q of the cantilever near 
the Rabi frequency Qr. Here, the quality factor of the un- 
loaded cantilever is Qc = 5 x lO"*, r=0.6 K. 



factor is gradually increased from about 200 to almost 
its bare value of 50,000 before it falls down to 5,000 at 
the Rabi frequency. The resonance frequency of the can- 
tilever (/c = 3 KHz) lies at the continuation of the right 
side of the curve shown in Fig. [71 The calculations show 
that at / = 3 KHz the quality factor Q w 1,800 with 
the frequency shift A/ w 0.32 Hz, which is in the range 
of the modified bandwidth (w 2 Hz). In principle, the 
dissipative part of the susceptibility Xzzi^) very sen- 
sitive to its parameters, especially, to 5, F^, and ili. If 
we had taken, for example, T2 = 1/F^ — 10~^ s, other 
parameters being unchanged, we would obtain Q « 200 
near 3 kHz, the resonance of the cantilever. Hence, while 
in the real experiment the external parameters, such as 
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FIG. 7: The modified quality factor Q of the cantilever near 
the Debye-like peak. Here, the quality factor Qc of the un- 
loaded cantilever is Qc = 5 x 10*, and T=0.6 K. 
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FIG. 8: Frequency shift, A^mXzz('^)/2wc, due to the back- 
action of the spins, versus the normalized frequency uj/^Ih, 
for T = 0.6 K. 

6 and fii can be controlled, the estimation of the damp- 
ing effect at the given frequency requires the knowledge, 
with good accuracy, of the spin-spin relaxation time T2 . 



B. Effect of the nuclear spin noise on the damping 
of the cantilever 

If the spin-lattice relaxation time Ti is extremely long, 
which happens at low temperatures, it is not possible to 
use conventional NMR methods, which rely on the mea- 
surement of the equilibrium spin polarization. In addi- 
tion, for nanoscale volumes (below about (lOOnm)^) the 
statistical spin polarization exceeds the mean Boltzmann 
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polarization^"'^. Hence, an alternative approach would be 
to measure a naturally-occurring statistical polarization: 
the spin noise^ii^i^. In this case, a sample has a mag- 
netic moment with a mean-squared- value proportional to 
^/N /X, where N is the number of nuclear spins in the res- 
onant slice, and /i is the magnetic moment of a single 
particleSai^. 

Hence, for long Ti we take for p{t'), p+{t'), and p-{t'), 
in Eqs. (|B8|, and (|BT0)) their initial values: p(0), 

p+{0) , and p-(0). These values correspond to natural 
spin fluctuations in a sample which is in thermal equilib- 
rium and under no external influence {coq = 0,Qi = 0). 
In order to find p(0), p+(0), and p-(0), we put Qi = 
in Eqs. (HH), and In this case the Rabi levels 

disappear and we get a sample in a constant polarizing 
field Bq with p^^^ = —TiT-, the equilibrium normalized 
population difference between two levels (see Eq. (|26p ). 
Hence, for this case, the quantity p*^"^ which before was 
the population difference between Rabi levels, remains 
the population difference between levels of a spin in the 
field Bq. Therefore, in the absence of an external field 
{Bq = 0) it is reasonable to consider p'"' as p{0), the ap- 
parent normalized population difference which provides 
the mean-squared-value of the naturally-occurring mag- 
netic moment p. Hence, 



1 



Np 



N 

Here, the quantity p(0) is negative because the upper 
level is now less populated than the lower one. 

In the same limit (fii = 0) we obtain from and 

(p9)) p^'' = p^-* = 0. This result does not depend on Bq 

and remains unchanged when Bq = 0. The quantities p^^"^ 

and p^'' describe the density matrix elements between 
Rabi levels (see Eqs. (26) in"*^) which should tend to 
zero when Rabi levels disappear. Hence, in the absence 



of the polarizing field Bq, it is reasonable to consider 
p(_") = p+(0) = and pL"^ = p_(0) - 0. 

Therefore, for the corresponding susceptibilities, we 
obtain the expressions (jBlSp . (|BI6p . and (|BI7p . Hence, 
in this case the expression for Xzz{i^) becomes: 



(46) 

where d{uj) is given in (|BI8p . 

The use of the spin noise for the MRFM 2D and 3D 
image reconstruction of the nuclear spin density was de- 
scribed in Refs.[l3 and0. These papers briefly reported 
that they measured an unexplained substantial decrease 
of the quality factor of the cantilever (from 50,000 to 
8,000 in Ref. [l^, and from 30,000 to several thousand 
in Ref. 0). Here we show that, qualitatively, this effect 
might be explained by the back-action of spin noise on the 
cantilever motion. 

For the quantitative estimate of this effect, we take as 
a guide the necessary parameters from Ref. [ij. Here we 
now assume that the spin-lattice relaxation time Ti is so 
long that it prevents the manipulation of the equilibrium 
Boltzmann polarization. For example, for the atom ^^F 
in calcium fluoride CaF2 studied in Ref. |T2, the time Ti 
at the experimental^ temperature T = 0.6 K was about 
10'^ s. 

The investigation of the imaginary part of the mag- 
netic susceptibility shows that its dependence on 
the high-frequency detuning S shows a sharp peak at the 
point of resonance, (5 = (see Fig. [TUl where the func- 
tion Xzzi^) is drawn at the resonance of the cantilever). 
Thus, we now write down the real and imaginary parts of 
the magnetic susceptibility for the exact resonance 
((5 = 0), where we expect the maximum effect: 



hVN\d{Lu)\ 



(47) 
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FIG. 11: The dependence of the dissipative part of the spin 
susceptibility on the distance x from the point of exact reso- 
nance, and here ijjc/2n — 3 kHz. 



5 X 10'^. Hence, for zero high-frequency detuning {5 = 
0) we obtain from Eq. the modified quahty factor 
Q « 0.3. This enormous reduction of the quality factor 
is primarily due to the large number of nuclei TV which are 
simultaneously at the exact magnetic resonance ((5 = 0). 

However, we stress that, strictly speaking, our esti- 
mates above cannot be considered as the only possible ex- 
planation of the cantilever damping observed in Refs. [T^ 
and[l. One of the reasons for this is that the long side of 
the cantilever in Refs. [H andd was perpendicular to the 
sample surface, which is different from the design shown 
in Fig. [T] For their design, the damping of the cantilever 
due to the spin-noise back-action is more sensitive (com- 
pared to our case) to the density distribution of the spins 
over the sample surface. Another reason is that here we 
assume that all resonant spins feel the same field and are 
located at the same distance from the tip just beneath it. 
Hence, to obtain more realistic values of Q, it is neces- 
sary to modify Eq. ([1]) with a more careful account of the 
density distribution of the resonant spins over a sample. 



and here 

\d{u)\' = co^ico^ -nl- Tir + rlinl 2coY ■ m 

From Eq. (l48l) it follows that the contribution of the spins 
to the damping of the cantilever is always positive, which 
means that the spin noise is cooling the cantilever. 

The width of the distribution A6 shown in Fig. [TO] is 
directly connected to the thickness of the resonant slice 
Ax = AS/'^G, where G is the magnetic field gradient. 
The dependence of the dissipative part of the spin suscep- 
tibility on the distance from the point of exact resonance 
is shown in Fig.[Tl]for = 3 T and G = 1.5 x 10^ T/m. 
From this figure, we estimate the effective thickness of 
the resonant slice as the full-width at half-maximum of 
the curve shown in Fig. 1111 Thus, we obtain Ax « 1 nm, 
which corresponds to the width of high-frequency detun- 
ing AS « 5 X 10~^wo (The full width at half maximum 
of the curve in Fig. [TO)) . 

For example, if we take for the sample the dimensions 
90 nm X 90 nm x 80 nm, with roughly 30 million nuclear 
spins (see Ref. [T2h . we obtain the effective number of 
spins in the resonant slice 90 nm x 90 nm x 1 nm that 
gives the main contribution to the effect: Nes ~ 3.5 x 10^. 

Below we estimate, from Eq. P5[) . the contribution 
of the spin noise to the cantilever damping at the res- 
onance frequency of the cantilever ujc, assuming T^p 3> 
u;/2n,ni/2'K. The result is (see Eq.[TH|: 

In order to estimate the modified quality factor from ([50]l . 
we take N w 10^, the number of spins in the resonant 
slice. Ftp w 10^ s^"'^, fc = 3 kHz, and the bare Qc = 



V. CONCLUSION 

In this paper we investigate the interaction of the 
MRFM cantilever with a system of nuclear spins. We 
show that the back-action of nuclear spins results in an 
additional contribution to the damping and frequency 
shift of the cantilever vibrations. We also show that a 
spin system may significantly change the quality factor 
of the cantilever. The cantilever can be either heated 
or cooled, depending on the sign of the high-frequency 
detuning. This effect exhibits a resonant nature, with 
a maximum at the Rabi frequency. We show that the 
main reason for this effect is that the longitudinal mag- 
netization (which is commonly measured in MRFM ex- 
periments) exhibits a resonance at the Rabi frequency, 
which can significantly alter its low-frequency evolution. 
We also analyze the influence of the spin noise on the 
cantilever damping and show that the spin noise may 
lead to a significant reduction of the quality factor of the 
cantilever. 

The interesting question is the lower bound on the cool- 
ing of the cantilever by nuclear spins. It might seem at 
first that the lower limit on the cooling of the cantilever 
is set by its zero-point fluctuations. However, this is not 
the case. From general considerations, the lower limit is 
set by the direct contribution of the spin noise to the can- 
tilever fluctuations (see, for example, Ref.l25l). In order 
to make a reasonable estimate of this limiting tempera- 
ture, it is necessary to: a) first calculate the spectrum of 
the spin fluctuations under a high-frequency field, and b) 
afterwards to consider a small number of Rabi photons, 
which requires treating the cantilever quantum mechan- 
ically. In our paper, we treat the spin system quantum 
mechanically, while the number of Rabi photons is large, 
which means that the cantilever behaves classically. This 
problem will be the subject of future investigations. 
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APPENDIX A: DYNAMICS OF AN IRRADIATED 
SPIN IN THE DRESSED STATE APPROACH 

Here we very briefly summarize some results of Ref . |43| . 
as applied to the problem we study in this paper. 

In addition to Eqs. (HH), and ([21]) (which describe 
the transition between Rabi levels of an irradiated spin) 
the (spin-|-field) system can also be characterized by the 
density matrix elements k's which describe the transi- 
tions between levels whith photon numbers that differ 
by one. These levels are approximately separated by hujQ, 
the energy between the levels of a bare spin. The rate 
equations for the k's are^: 



These quantities are directly connected to the longitudi- 
nal magnetization of a sample with N spin- 1/2 particles: 

Mz = -^{<Tz), (A8) 



and its transverse components 

Mx = ^{<^x) ; My = ^ (ay) 



(A9) 



It is very instructive here to show that the steady- 
state solution for the density matrix provides the well- 
known Bloch expressions for the longitudinal and trans- 
verse components of the magnetization. 

(st) 

The stationary magnetization , which is defined 



in is obtained from ^ with the help of ((?7)) . 

By using the substitutions = I/T2, fli = jBi, we 



obtain for M 



(st). 



M^z^^ = Mo 



1 + {T26Y 



l + iT2df + TiT2{-fBif 



(AlO) 



where 



Mo^^tanhfMi 
" 2 \2kBT 

Now we find the steady-state solutions of Eqs. (|Alll . (|A2|1 . 
^31) , and (IA4D . It is not difficult to see that the solution 



of these equations has the form: k+ = c+e 
ce~'""™*, K+ — c+e"*""™*, K_ ~ c_e^''^'"™*, where 



c = c+— +r!2„)cos20 



(All) 



It 



dK+ 



(Al) 



-iuj^ujK - AiK + Bk+ + K+{T_) cos2d , (A2) 



c_ = — zc^-^ri/jFip sin 20 , 



(A12) 
(A13) 



~dr 



-iujmwi^+ — i^Ri^- +_Bk — + K (r_) sin 29 , 

(A3) 



-iuJmwH- — i^RH+ ~ ^ipH- ■ (-A-4) 



dt 



For possible applications of this method, the quantities to 
be measured are the averages of the Pau li sp in operators 
{<7x), (cy), (cTz). As was shown in Ref. l43l 

(az) = p{t) cos 29 + p+ (t) sin 29 , ( A5) 



(ax) = sin2e'Re[K(t)] - cos 26'Re[K+(i)] - Re[K_(i)], 

(A6) 

(cry) = -sin26'Im[K(t)] + cos 26'Im[K+(<)] +Im[K_(t)], 

(A7) 



(A14) 



D = 

Ti 

By choosing c+ = 1 and using Eqs. (|A6|) and (|A7|) to 
calculate (IA9I) we find 



Mx — Mc costJniwi + Ms sinwmwi , (A-15) 
My = —Ms cos Wmw^ + Mc sin uo^^jt , (A-16) 



where 



Me = Mo 



Ms = Mo 



l + iT2Sf +TiT2{-fBif 



T21B1 



l + {T25f +TiT2{-iBif 



(A17) 



(A18) 



Therefore, the steady-state solution of the density ma- 
trix equations (HHl-lllTI) and (|M|) - (|X4)) provides the weU 
known expressions for the longitudinal (jAlOp and trans- 
verse components (jAlTI) , (|A18I) of the magnetization^. 
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APPENDIX B: LOW-FREQUENCY RESPONSE ( ^ A-^ ^P+i^) _ ^ _^ ■^^^ '^P-ft) 



OF AN IRRADIATED TWO-LEVEL SYSTEM \dt J Sf{t') 6f{t') Sf{t') 



The response of a two- level system to a weak external ~ h )P-i )cos , ( ) 

force f{t), in the limit of vanishing force f{t), is found 
from the following equations: 



^ = -Aip + Bp+ - yit)p-sm2e + ir-)cos29 , , \ x U\ A U\ 



^ - -i^RP- +Bp- A2P+ + y{t)p- cos 29 
at h 

+ (r_)sin26', (B2) 



-S(t - t')\p+ {t') cos 20 - p(t') sin 26*1 . (B6) 
n 



where we used the definition 5f{t)/5f{t') = 5{t — t'). 



The corresponding susceptibilities Xp{^)i Xp+(^)j 
7 . Xp-{^) ^re defined similar to Eq. ^TU\i : 

= -iniip+-r^p_ + ^f{t) {p+cos2e - psm29) . 

(B3) 

By taking the functional derivative with respect to /(f') Sf{t') ^ J 27r ^"^^ ^ ^ Xp('^)j 6tc. (B7) 

we obtain 

Vdt / '^/l* i '^/l* ) From Eqs. jll}, HU- and JB6| wc can readily find 

(B4) 

the susceptibilities Xp(^)i Xp+C"-^); and Xp_((jj): 



XpH--^^-(0sin2^ 



p+(t')cos20-p(t')sin26' } , (B8) 



Xp+ (^) = A p- cos 29 i~iuj + T^){-iLu + Ai)- p^ {t') sin 263 {-iuj + T^) 

p+it')coa2e ~ p{t')sm29 {-iuj + Ai) | , 



(B9) 



Xp- (^) 



hD{Lu) 



(^-iu;+ ^^{-iu + r^) p+{t') cos 29 - p{t')sin2e - cos26'| , (BIO) 



where 

D{uj) = {-iuj + r^)^ ^-iw +y)+ {-i^ + ^1) f^ls • 

(Bll) 

The functional derivatives (|B7|) are defined for t > t', 
where t' is the time the external force is being applied. 
Hence, the susceptibilities (|B8p . (jB9|, and (jB10[) describe 
the evolution of the system for times t > t' , with the 
p{t'Ys, corresponding to when the external force is ap- 
plied. Therefore, the subsequent evolution of the system 



depends on its state just before the perturbation is ap- 
plied. 

In what follows we consider two cases. The first one is 
when the relaxation time Ti is relatively short. In this 
case the system quickly reaches thermal equilibrium dur- 
ing the measurement. For this case, we take for p{t'), 
p+{t'), and p-{t'), their steady-state values: from 
Eq. (EH), pff*^ from Eq. and from Eq. For 

the corresponding susceptibilities, we obtain the follow- 
ing expressions^: 
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02 

-^B cos 29 



(B12) 



/ \ on (0) 



/ 1 ^ ^1 

{-iuj + r^) ( -iu + — ) - {-i^ + ^1) jr- 



Ti 



(B13) 



hDiLo)ri 



(B14) 



The other case is when the spin-lattice relaxation time Ti 
is extremely long compared to the measurement time. If 
we measure the spin noise in this case, then it is reason- 
able to take p(0) = -1/VN, p(0)+ = p{0)_ = (see the 
explanation in Section III). For the corresponding sus- 
ceptibilities, we obtain from (|B8l) . (|B9)) . and (jBlOp (in 
the limit 1/Ti < uj): 



where 

d(u;) = {nj - 2^2) + (tj2 _ f72^ _ ) . (B18) 
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